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As the geosynchronous orbit for satellite communication becomes 
increasingly crowded, the effect of adjacent satellite interference is of 
primary concern. We describe here the saddle-point approximation 
technique for evaluating the error probability of M-ary PSK systems 
with adjacent satellite interference. In comparison with previous 
methods, this technique has the advantages of both computational 
simplicity and accuracy. Results include a sample calculation to 
study the effect of fading on the uplink of a 12/14-GHz quadraphase- 
shi ft- keying system with 3-degree satellite spacings and 3-m earth 
stations. 

I. INTRODUCTION 

The problem of adjacent satellite interference, which is basically one 
of co-channel interference, has received extensive attention in terms 
of exact computation methods 1 " 5 and bounds. 6 " 9 In spite of this, system 
designers still often treat the interference as additive Gaussian noise, 
even though this has been shown to be very pessimistic 2,4 for lower 
error rates (<10" 8 ). The reason for this is that the exact methods are 
too difficult to use and the simple bounds are not tight enough for 
large signal-to-noise ratios (i.e., low error rates). As a result, the 
communication system is overspecified and the geosynchronous orbit, 
a limited natural resource, is not used efficiently. 

In this paper, we describe a saddle-point method 1011 which is both 
simple to use and accurate for evaluating adjacent satellite interfer- 
ence. After defining the problem in Section II, we derive the saddle- 
point approximation from the Fourier inversion integral representation 
of the error probability in Section III. To illustrate the accuracy of the 
approximation, we consider the uplink error performance of a QPSK 
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system in the 12/14 GHz band, with 3m earth stations and 10 interfer- 
ing satellites spaced 3 degrees apart. 

II. PROBLEM STATEMENT 

We consider a M-ary PSK system with coherent detection. The 
transmitted signal is one of M possible waveforms: 

s m (t) = y/2S cos(wo* + 6 m ), < * < T, (1) 

where 6 m = 7r(2m — l)/M is the carrier phase angle corresponding to 
the mth symbol (m = 1, 2, •••, M). At the receiver, there are K 
interfering signals from different PSK transmitters: 

K 

I(t) = E^ cos(w,* + to + <f>,), (2) 

where /, and uj are constants, »//, is a random variable conveying the 
information from the y'th interferer (j = 1, 2, ■ ■ • , K), and tfy is a 
uniformly distributed random variable in [0, 2tt). Then for s m (t) trans- 
mitted, the received signal is given by 

r(t) = s m (t) + I(t) + n(t), (3) 

where the receiver noise, n(t), is a stationary, zero-mean Gaussian 
process with variance a 2 . Rewriting (3) in quadrature components, we 
obtain 

K 

y/2S cos m + X Jzij cos \j(t) + n c {t) 



r(t) = 



/-J 



COS O)ot 



\/2S sin m + X -J2Ij sin Xj(t) + n s (t) 



sin a>ot. 



The receiver samples the envelopes of the in-phase and quadrature 
terms at t = to, and chooses one of M symbols as the transmitted one. 
Without loss of generality, we assume 6\ was sent. The envelope 
samples are given by 

K 

r c = ^/2S cos 6i+ J Jzij cos A, + n c (4) 

and 

K 

r s = y/2S sin ft + £ ^27/ sin Ay + n s , (5) 

where the A/s are independent and uniformly distributed in [0, 2tt), 
the n c and n a are zero-mean Gaussian with variance a 2 , and the 
dependence on to has been suppressed. If the phase angle of the vector 
(r c , r s ) lies within [0, 2tt/M), the receiver would decide that 0\ was 
transmitted. As shown in Ref. 4, the symbol error probability is 
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bounded by 

P e < 2Pr[r„ < 0]. (6) 

III. THE SADDLE-POINT APPROXIMATION 

The probability P s 4 Pr[r s < 0] can be evaluated via the Fourier 
inversion integral: 12 

Ps = -^-. I Q(u)du/u, (7) 

2m Jc 

where the path of integration C is along the real axis (— oo to oo) except 
for an indentation above the origin, and O(m) is the characteristic 
function of r s / V2S: 

K 



with 



and 



O(h) = exp[iW - u 2 /4p 2 ] I] MRju), 



d= sin(ir/M), 

p 2 = S/a 2 , the carrier-to-noise ratio (cnr), 
Rj = yfljJS, they'th interference-to-carrier ratio, 

Jo{x), the zeroth order Bessel function of the first kind. 



Define G(u) = ln[0(«)]. Then the saddle points of the integrand in 
(7) (ignoring the effect of u~ l ) are solutions of G'(u) = 0, i.e., 



K 



id - u/2 P 2 - I RjJi(Rju)/MRju) = 0. (8) 

In Appendix A, we show a unique solution of (8) on the positive 
imaginary axis, i.e., there is only one saddle point. Denoting this by u s 
= iy s , we can rewrite the saddle-point equation and define v s as the 
solution to 

H (y) = -d + y/2 P 2 + I RjIi(Rjy)/URjy) = 0, (9) 

y=i 
where I (x) and I\{x) are the zeroth and first-order modified Bessel 
functions, respectively. Since it is observed in Appendix A that 
H{0)H"(0) > and H'( y) and H"(y) do not change sign in [0, y s ], 
Newton's root-finding method (Ref. 13, p. 18): 

y k+i =y k -H(y k )/H'(y k ) (10) 

would converge monotonically with an initial estimate of ya = 0. 
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Since 3>(a) is analytic in the entire complex u-plane, we can move 
the contour of integration for (7) away from the pole at the origin, to 
intersect the imaginary axis at u s . Then expanding G(u) around the 
saddle point, we obtain 

G(u) = G(u s ) + V 2 (u - u s ) 2 G"{u s ) + • • • 

A G 8 -%(« - u s ) 2 a 2 s + ■•-. 
Substituting this Taylor series expansion of G(u) into (7), we can 
integrate term by term to obtain an asymptotic series (as is done in 
Appendix B), with the leading term given by 

P s = V 2 exp[G s + 3«/2]erfc[WV2]. (U) 

where erfc(^) is the complementary error function (Ref. 13, p. 297). 
We shall refer to P s as the saddle-point approximation of P s . It is 
shown in Appendix B that, as y s — > oo, which corresponds to either 
P 2 (cnr) — *■ oo for an open-eye pattern (d > D =£>-i Rj), or D — d — » 
for a closed one (d < D), the saddle-point approximation is asymp- 
totically exact. More precisely, 

P S = P S + IK exp[G s ]/ y/frr 8a 5 s y 5 s + 0(y7 6 ). (12) 

Fortran subroutines for determining the saddle point and evaluating 
P s are given in Appendix C. 

The deformation of the path of integration to pass through the 
saddle point provides not only an accurate approximation, but also a 
very efficient method of evaluating the inversion integral numerically. 14 
Let u = x + iy 8 . Then eq. (7) becomes 



/ 



P a = exp[G s ] ) exp[-x74p z + ix(d - y s /2p 2 )] 

K 

• II {Jo[Rj(x + iy s )]/Io[Rjy s ]) dx/27r(y s - ix). (13) 
y=i 

Using the inequality 

\Jo(x + iy)\<I (y), (14) 

we immediately obtain a simple upper bound in terms of .y s : 

P s < exp[G s ]p/y s v£ A P u , (15) 

which is tighter than the Chernoff bound 15 for .y s > p/V^r. To truncate 
the range of integration in (13) to (—X, X), we first bound the remain- 
der: 

\R X \ < exp[G„] exp[-x 2 /4p 2 ]dx/7ry s = P u erfc[X/2p] 
Jx 

< P u exp[-X 2 /4p 2 ]. (16) 

1 142 THE BELL SYSTEM TECHNICAL JOURNAL, SEPTEMBER 1 980 



Then the truncated integral is given by 

P B - R x = exp[GJ exp[-* 2 /4p 2 + E r (x)] 

Jo 

• { v s cos Ei(x) - x sin Ei(x)} dx/ir(x 2 + yl), (17) 
where 
exp[£ r (x) + iEi(x)] 

K 

= exp[ix(d - y a /2 P 2 )] U {J [Rj(x + iy s )]/Io[Rjys]}. 

Therefore, with (16) and (17), P a can be determined to any degree of 
accuracy desired, and (17) is easier to compute numerically than (7) 
because the latter has a singularity at the origin. 14 

IV. EXAMPLE 

To illustrate the accuracy of the saddle-point approximation, we 
consider the case of a geostationary communication satellite, flanked 
on each side by K/2 interfering satellites equally spaced by A0 degrees. 
According to the model in Ref. 4, if all satellites transmit at the same 
power level, the yth interference-to-carrier ratio at the earth station is 
given by 

Rj r- (kM)~ 1M 18\/D a , where k = \ *■ {zL^L (18) 

[ Ji] > A/Z 

D a /X is the ratio of diameter to wavelength of the earth station 
antenna, and is approximately 100 for 3m stations in the 12/14 GHz 
band. 

If K + 1 earth stations are transmitting to their respective satellites 
with the same ccir antenna sidelobe characteristic, 4 the interference- 
to-carrier ratios at the desired satellite are also given by (18) (see Fig. 
3). We assume that the earth stations are separated far apart geograph- 
ically to experience different fading conditions due to precipitation, 
but close enough so that the angle subtended by the maximum sepa- 
ration at the satellite is small as compared to A0. Then the worst case 
condition is when the desired uplink signal is subject to fading while 
the interfering signals are not, and (18) is modified as 

Rj = qRj = 4>(kdO)- V3S 18X/D a , (19) 

with the fading parameter, 1 < <#> < ». 

Figure 1 shows the uplink performance of a QPSK (M = 4) system 
with interfering satellites at A0 = 3° spacings and 3m earth stations, 
for 0-, 6-, and 12-dB fading. The two parameters of interest are the 
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10-3 



10- 



10~ 5 — 



10" 6 



12dB FADE 
CIR = 10.7dB 
D/d = 1.02 



6dB FADE 

• CIR = 16.7dB 

D/d = 0.51 



OdB FADE 

CIR = 22.7dB- 

D/d = 0.26 



EXACT 



SADDLE POINT 

• UPPER BOUND 



20 



25 



CNR IN DECIBELS 



Fig. 1 — Performance of a QPSK system subject to fading, with 10 interfering satellites 
at 3-degree spacings. All curves are upper bounds on the actual symbol error rate [eq. 
(6)]. The earth station antenna has diameter 3m (D a = 100A). cut is the total carrier-to- 
interference ratio and D/d is the peak distortion. 



total carrier-to-interference ratio (cir): 

/ « IK/2 

p? = 1 / 2 Rf = (A9) 2 HD a /\18<t>) 2 / I 2k~ 25 , 

I y=i / k=i 

and the peak distortion ratio 



(20) 



We plot the upper bound (6) on the symbol error probability using P 8 
and P u from (11) and (15), respectively, and the exact P s using (16) 
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and (17). For fading conditions of and 6 dB, the exact and the 
approximate curves are indistinguishable. In the 12-dB case, we note 
that D/d > 1, which means the eye pattern is closed and the saddle 
point location does not approach infinity for large cnr. However, even 
though the saddle-point approximation is not exact asymptotically, it 
is still quite accurate, deviating from the exact result only in the second 
significant figure. 

To determine the fade margin against interference at the satellite, 
we plot cnr vs cir for different error rates in Fig. 2. For example, if 



SADDLE POINT 



GAUSSIAN 




= 10-6 



0-5 



20 

CIR IN DECIBELS 



30 



Fig. 2— Effect of variations in cir on the required cnr to maintain constant error 
rate. The results are applicable to the uplink QPSK system with 10 interfering satellites 
equally spaced on both sides of the desired one, regardless of orbital spacing and earth 
station antenna size. 
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^np <A? C ^P > < ^^ 

E_, E E, E 2 

Fig. 3— Equally spaced satellites and their earth stations. Earth station Ek (k j* 0) is 
pointed at satellite S*, but its signal would interfere through its sidelobe with the 
transmission from Ea to S . 



the system is designed for P e = 10~ 6 and a nominal cm of 23 dB, a 6- 
dB fade would require an additional increase of 1.7 dB in cnr, which 
means a total fade margin of 7.7 dB. For comparison, we also show the 
result from the Gaussian approximation in Fig. 2, using 

P s = x h erfc(dpe), 

with the equivalent cnr given by p7 2 = p~ 2 + pi 2 . The required fade 
margin in this case is 8.4 dB, which is too high by 0.7 dB; for lower 
nominal cir's, the discrepancy is even larger. Moreover, the Gaussian 
approximation would overestimate the niinimum cir for P e = 10~ 6 by 
at least 3.5 dB. 

As seen from (20), given a total cir, the allocation of interference 
power among the 10 interferers is independent of satellite spacing 
(A0) and the size of the earth station antenna (D a /\), so that Fig. 2 is 
applicable to the model in Ref. 4 for all satellite spacings and antenna 
sizes. 
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V. CONCLUSION 

In this paper we have described the saddle-point approximation for 
analyzing M-ary PSK systems with co-channel interference and dem- 
onstrated its usage and accuracy in an example. Along with the results 
for intersymbol interference," this method can be extended to analyze 
the combined effect of filtering, co-channel and adjacent channel 
interference. 16-18 It should be noted that the inequality (14) is a 
fundamental property of characteristic functions and that the upper 
bound (15) in terms of the saddle point is applicable to all types of 
additive interference, as long as the receiver noise is Gaussian. 
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APPENDIX A 

Uniqueness of the Imaginary-Axis Saddle Point 

The saddle-point equation (9) is repeated here for easy reference: 

H{y) = -d + y/2p 2 + £ RMRrf/URjy) = 0. (21) 

Define B(x) = Ii(x)/I Q (x). Then the derivative of H{y) is given by 
H'(y)=-X+Y t RfB'(R J y). 

Z P 7=1 

Since H(0) = —d and //(oo) = oo for p < oo, at least one solution to (21) 
exists in < y < oo. It can be shown that B'(x) > for < x < oo, so 
that //' ( y) > in < y < oo, which means that only one saddle point 
is on the positive imaginary axis. 

Figure 4 shows B(x) and its first two derivatives. We note that 
B"(x), and therefore H"(y), is strictly negative for positive arguments. 
Using the asymptotic expansion of I (x) and Ii(x), we obtain for large 
x, 

B{x) ~ 1-jc _ V2+ -.-, 



and 



B'(x) ~x-72 + *-74+ ••-, 

B"(x) ~-x- 3 -3x- 4 /4+ ••-. 
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Fig. 4— Behavior of B(x) = I,{x)/I„{.x) and its derivatives. 



APPENDIX B 

Asymptotic Expansion of P s 

Along the path of integration which is parallel to the real axis, 
passing through the saddle point u s = iy s , G{u) can be written as 

G(x + iys) = Fiy* - ix) 

= F(y a ) + (-ix)F'(ys) +i^V(*) + l ^V>(y s ) 



2! 



k\ 



AGs' — ol + FM + iFiix), 



with 



.2* 



F r (x)- I (-) l^P m M 
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and 



Fi(x) = I <-)' 



'(2A-1 



(*). 



A (2A-1)! 

Substituting (22) into the inversion integral (7), we obtain 



(22) 



P a = exp[G s ] exp 



— 61 + FM 



exp[iFi(x)]dx/{y s — ix)2ir 



= exp[G 8 ]y s exp 



-—0 2 s +F r (x) 



• \ cos Fi(x) - — sin F t (x) \ dx/(x 2 + y 2 )ir 



A exp[G s ]j s exp 



-T* 



1+ £ C k x 2k \dx/(x 2 +yi)7r, 
*=2 ) (23) 



with the first two C*'s given by 

C 2 = F (4) ( yil )/4\-F i3) (ys)/y s 3\ 
and 

C 3 = F {5) (y s )/y s 5\ - F {6) ( v s )/6! - (F (3) /3!) 2 /2. 

Then, with t = x 2 o 2 /2 and z = o 2 y 2 /2, term-by-term integration of 
(23) gives 



P s = Viexp[G s ] e~' 



i + y c*(2/ff5)V 



t~ 1/2 dt/(t + z)2tt 



= i exp[G s + z] j erfc(Vi) + - £ C*(2/aft* 

2 ^ 7T *=2 



■rlA + jWQ-*.*)}. 



(24) 



where T(a, 2) is the incomplete gamma function (Ref. 13, p. 260). 
To illustrate the asymptotic behavior of P s , for large y s , we have 

F"'(y) = -K/y 3 + 0(y- 4 ), 



and 



so that 



F {4) (y) = 3K/y 4 + 0(y" 5 ), 

F m (y) = 0(y- k ), for k > 5, 

C 2 = 7K»24! + 0( jj 5 ) 
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and 

C k = 0(y7 2k ), for k > 3. 
Then the discrepancy between P s and P a is given by 

Ps~Ps= yi * —Z 2 + 0(2 3 ). 

2Vtt-2 32 

However, it should be noted that y x may not approach infinity as 
p 2 — » co. For large y s , the saddle-point equation (9) becomes 

H(y a ) ~ y s /2p 2 + D- K/2y a - d = 0, 

where Z) = £/L] ^>» an ^ we have two asymptotic solutions for^ s : 

y s ~ 2p 2 (rf -D)+ K/2(d - D), for d > D 

and 

y s ~K/2(D-d), for rf < Z). 

Therefore j s would approach infinity either as p 2 becomes large in the 
open eye-pattern case (d > D), or as D - d tends to zero; if the eye 
pattern is closed {d<D), y s would approach a finite limit. 

APPENDIX C 
Fortran Subroutines 

The saddle-point approximation is computed to a desired accuracy 
using two subroutines. The number of interferers K, the individual 
interference-to-carrier ratios {/?(/)}, and d = sin(-n/M) are specified 
via the common statement. Subroutine "ggg" is used to calculate 

K 

g0 A G(iy) = -yd + y 2 /4 P 2 + I ln[I n (R jy )l 

K 

g\ A -iG'(iy) =H(y) = -d + y/2p 2 + % RMRjy)/h(Rjy), 

and 

g2 A -G"(iy) = H'{y) 

= ^p 2 +Z Rf{l\(Rjy)io(Rjy) - ftifl&Wttftjy)- 

The modified Bessel functions My) and I\(y) are evaluated by calling 
besri, and I\{y) = I (y) - I\(y)/y. 

Subroutine "sadpnt" implements the saddle-point search algorithm 
as defined in (10), terminating when the latest relative change in the 
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saddle-point estimate is less than eps. It then computes P s according 
to (11). 

subroutine sadpnt(ps, eps) 

common/ ggg/y, d, ro2, K, R(20) 
c ps — saddle point approximation 
c eps — relative accuracy of saddle point estimate 
c All other parameters are transmitted via the labelled common. 

y = 0.0 linitial guess 

5 caQ.ggg(g$,gl,g2) 

7 yinc = -gl/g2 

y = y + yinc 

if(abs(yinc/y) > eps) go to 5 

caUggg(g0,gl,g2) 

ysig = sqrt(fl.5 * g2) * y 

ps = 0.5 * exp(g0) * erfw(ysig) !erfw(2) = exp(z*z) * erfc(z) 

return 

end 



subroutine ggg(g$, gl, g2) 
c Constant amplitude sinusoidal interferers 
common/ ggg/y, d, ro2, K, R (20) 
dimension 6(2) 
sig2 = .5/ro2 

gfy = 0.5 * v * v * sig2 — v * d 
gl = y * sig2 — d 
g2 = sig2 
do 10y= 1,K 
rj = R(j) 
ry = rj*y 
nb = 2 
call besri(ry» rib, b) IModified Bessel fens: Ik, k < nb 

60 = 6(1) 

61 = 6(2) 
blp = 0.5 

if(ry.ne.0.) 61p = bl/ry 
blp = 60 - blp 

g$ = g</> + alog(bfi) 
gl=gl + rj * 61/60 

g2=g2 + rj * rj * (blp * 60 - 61 * 61)/(60 * 60) 
10 continue 

return 
end 
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